Search for: All records

Many researchers have used the birefringence of P‑to‑S converted waves from the Moho discontinuity to constrain the anisotropy of Earth’s crust. However, this practice ignores the substantial influence that anisotropy has on the initial amplitude of the converted wave, which adds to the splitting acquired during its propagation from Moho to the seismometer. We find that large variations in Ps birefringence estimates with back-azimuth occur theoretically in the presence of P‑wave anisotropy, which normally accompanies S‑wave anisotropy. The variations are largest for crustal anisotropy with a tilted axis of symmetry, a geometry that is often neglected in birefringence interpretations, but is commonly found in Earth’s crust. We simulated globally-distributed P‑coda datasets for 36 distinct 4‑layer crustal models with combinations of elliptical shear anisotropy or compressional anisotropy, and also incorporated the higher-order anisotropic Backus parameter C. We tested both horizontal and tilted symmetry-axis geometries and tested the birefringence tradeoff associated with Ps converted phases at the top and bottom of a thin high‑ or low‑velocity basal layer. We computed composite receiver functions (RFs) with harmonic regression over back azimuth, using multipletaper correlation with moveout corrections for the epicentral distances of 471 events, to simulate a realistic data set. We estimate Ps birefringence from the radial and transverse RFs, a strategy that is similar to previous studies. We find that Ps splitting can be a useful indicator of bulk crustal anisotropy only under restricted circumstance, either in media with no compressional anisotropy, or if the symmetry axis is horizontal throughout. In other, more-realistic cases, the inferred fast polarization of Ps birefringence estimated from synthetic RFs tends either to drift with back-azimuth, form weak penalty-function minima, or return splitting times that depend on the thickness of an anisotropic layer, rather than the birefringence accumulated within it.   

Abstract Abrupt velocity gradients in the upper mantle, detectable by receiver functions (RF) techniques, have been known to exist down to the depths of ∼110 km beneath northeastern North America. Comparisons with the surface wave velocity models have designated some negative velocity gradients (NVGs) as the lithosphere‐asthenosphere boundary (LAB), delineating a relatively thin lithosphere beneath this region. This work presents a systematic survey of upper mantle layering in seismic properties using P‐S RF analysis at 62 long‐running sites with dense lateral sampling. We examine both radial and transverse component RF for indicators of seismic anisotropy and adopt the notion of seismic attributes, utilized in active‐source seismology, to characterize the spatial distribution of directionally variant and invariant signal components. We confirm a widespread presence of NVGs at depths 60–100 km throughout the region, consistent with previous studies using mode‐converted body waves. We also find numerous converting boundaries that reflect changes in directional variation (anisotropy) of seismic velocity, indicating complexity of rock texture in the upper mantle. Some of these boundaries appear as deep as 185 km, implying that the lithosphere extends much deeper than the widespread NVGs would suggest. In this, our results agree with recent estimates of the lithospheric thickness in thermodynamically consistent models combining seismic, gravity, and heat flow constraints. 

Abstract We use splitting in core‐refracted teleseismic shear waves (SKS, PKS, and similar) to investigate anisotropic properties of the upper mantle beneath the Superior craton in eastern North America and the Yilgarn craton in Western Australia. At four sites in each craton, we assemble extensive data sets that emphasize directional coverage, and use three different measurement methods to develop mutually consistent constraints on the nature of splitting and on the likely anisotropic properties that cause it. In both cratons, we see evidence of clear directional variation in both delays and fast polarization directions, as well as lateral differences between sites. Relatively small (0.3–0.8 s) amounts of splitting imply weak anisotropy within 150–220 km thick mantle lithosphere. Anisotropy in the asthenosphere likely contributes to splitting in North America where fast directions align with absolute plate motion, but not in Western Australia where fast polarizations and plate motion are nearly orthogonal. 

Abstract Layering within the cratonic lithosphere has been explored and reported in different cratons using a range of techniques. However, whether there exists a common feature in the lithosphere for all the cratons is not clear yet. In this study, we carry out a comparison study between the Yilgarn craton in Western Australia and the Superior craton in North America that have never been in direct contact throughout their tectonic history. To have a detailed description of the lithospheric layering in both cratons, we employ receiver function analysis with harmonic decomposition to characterize the anisotropic seismic structure beneath 4 long‐operating sites in each craton. We can identify multiple unique anisotropic boundaries above 170 km at all sites in both cratons. Properties of the anisotropic boundaries are distinct both within and across the cratons. Our observation agrees with a commonly accepted view of the cratonic lithosphere consisting of at least two layers. Moreover, it adds new details to the previous view and reveals lateral variations of the anisotropic properties over distances of a few hundreds of kilometers. Such variations in anisotropic properties likely reflect the tectonic history predating the final assembly of cratons, and suggest horizontal movements are necessary for the formation of cratonic lithosphere.